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ABSTRACT 



We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to 

mm-wave photometry of quasars with redshifts < z < 6 and absolute /-band magnitudes -22 > /ab > -32 We use the science demonstration 

observations of the first ~ 16deg^ from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars 

^>^ _ from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data 

\^ , with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area 

at > 5a- at 250,350 and 500 //m. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: 

low-luminosity quasars with absolute magnitudes in the range -22 > /ab > -24 have a comoving star formation rate (derived from 100 /jm 

rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with /ab < -26 have a maximum contribution to 

Q , the star formation density at z ~ 3. The volume-averaged star formation rate of -22 > /ab > -24 quasars evolves as (1 + i)^-^^oj gj ^ < 2, but 

$^ ■ the evolution at higher luminosities is much faster reaching (1 + z)'"*' at -26 > /ab > -28. We tentatively interpret this as a combination of a 

"^ ■ declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation. 
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'^ 1. Introduction star formation in their host galaxies (e.g. Granato et al. 2006), 

^-^ ' but these feedback models have many adjustable parameters. 

^ . The cosmic star formation history (e.g. Madau et al. 1996 and Feedback is arguably the principal uncertainty in semi-analytic 

^ ; others) was quickly realised to bear a striking apparent siim- ^^^^^^ ^f g^l^^y evolution. One of the few observational ap- 

';^. . larity to the evolving luimnosity density of quasars (QSOs) at p^oaches available to constrain feedback models is to measure 

in most redshifts (e.g. Boyle & Terlevich 1998, Franceschini et al. ^^e star formation rates in QSO host galaxies, but so far very 

O . 1999), suggesting a link between the physical drivers of black f^^ q^q^ ^^^^ ^^^^^^ far-infrared, submm or mm-wave detec- 

O . hole growth and stellar mass assembly. Other data also indirectly ^^^^^ ^^^^ ^^^^^ ^j^. formation rates could be inferred. 
'""' • suggested links. Mid-infrared spectra (e.g. Genzel et al. 1998, 

j>! [ Spoon et al. 2007) and radiative transfer modelling (e.g. Farrah The Herschel Astrophysical Terahertz Large Area Survey 

K~J et al. 2002) imply higher luminosity starbursts have increasingly (H-ATLAS, Bales et al. 2010) is the largest open time key 

/\ large bolometric fractions from active galactic nuclei (AGN). project on the Herschel Space Observatory (Pilbratt et al. 2010). 

H The tight /T-band Hubble diagram of far-infrared-selected hyper- The survey aims to map 550 deg^ at five wavelengths from 

. - - : luminous starbursts also bears a striking similarity to the K - z 110-500;umto5crdepthsin the range 32-50 mJy at > 250 ;um. 

relation of radiogalaxies (e.g. Serjeant et al. 2003, though see Among the key science goals is a constraint on the star forma- 

Smail et al. 2004 for their submm-selected counterparts). tion rates of many thousands of QSOs through far-infrared and 

A close relationship between black hole growth and stellar submm photometry. In preparation, Serjeant & Hatziminaoglou 

mass assembly is also demanded by the observed close correla- (2009) used a compilation of available far-infrared to mm-wave 

tions betwen supermassive black hole masses and spheroid prop- photometry of QSOs from IRAS, ISO, Spitzer, SCUBA and 

erties (e.g. Magorrian et al. 1998, Merritt & Ferrarese 2001), MAMBO, to predict the numbers of QSOs directly detectable by 

which exist despite the enormous disparities of spatial scales and H-ATLAS. Several hundred QSO direct detections are expected 

masses. Such correlations are predicted by models with radia- in H-ATLAS and the first H-ATLAS detections are described in 

tive and/or kinetic energy outputs from the AGN regulating the e.g. Gonzalez-Nuevo et al. (2010). Serjeant & Hatziminaoglou 

(2009) also assumed an M82 spectral energy distribution (SED), 

* Herschel is an ESA space observatory with science instruments used observations close to 100(1 + z) pm where possible, and 

provided by European-led Principal Investigator consortia with impor- found the stacked lOOyum luminosities of QSOs typically scal- 

tant participation from NASA. ing roughly as the square root of the optical luminosities with 
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a redshift-dependent normalisation, supporting the idea of cou- 
pled black hole mass and stellar mass assemblies. Almost iden- 
tical results were obtained with an Arp220 SED. This SED in- 
sensitivity can be readily understood: at a fixed bolometric lu- 
minosity, the SWIRE template SEDs (Polletta et al. 2007) for 
the starbursts M82, Arp 220, NGC 6090, IRAS 2055 1 -4250 and 
IRAS 2249 1 - 1 808 have 1 00 fj.m monochromatic luminosities all 
within a factor of two (note that we are only concerned with the 
starburst bolometric contribution). 

In this paper we extend this analysis to the first data from 
H-ATLAS. We interpret the lOOyum luminosities as star forma- 
tion, since AGN dust tori SEDs are expected to peak in the 
mid-infrared, (e.g. Efstathiou & Rowan-Robinson 1995, though 
see the discussion in Netzer et al. 2007). We estimate the mean 
star formation rate in bins of redshift and absolute /ab magni- 
tude, then use the QSO luminosity function to make the first 
constraints on the star formation rate per comoving volume of 
QSOs from < z < 6 and a factor of lO** in optical luminosity. 
We assume density parameters Qm = 0.3 and Oa = 0.7 and a 
Hubble constant of //q = 70kms~' Mpc '. /-band QSO magni- 
tudes are quoted in the AB system with /iT-corrections assuming 
dln5v/dlnv = -0.5 and no internal dust extinction correction. 
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Fig. 1. The distribution of QSOs considered in this paper in 
the optical luminosity-redshift plane. Red crosses mark QSOs 
in the Herschel ATLAS survey. Open circles are SDSS QSOs 
with SWIRE coverage, while filled circles are Palomar-Green 
QSOs with IRAS and B-band data. Diamonds represent QSOs 
observed at 850 /im while open squares represent QSOs ob- 
served at 1200/im. Adding H-ATLAS QSOs extends the range 
of absolute magnitudes at redshifts z < 3. 



2. Data acquisition and analysis 

2.1. Herschel ATLAS 

We use only data from the SPIRE instrument (Griflin et al. 2010) 
in this paper For more details of the data analysis see Pascale 
et al. 2010 in prep.; we summarise the main points here. The 
SPIRE images were registered to a common reference frame us- 
ing stacking analyses of Sloan Digital Sky Survey (SDSS) galax- 
ies. Neptune was used for flux calibration, requiring multiplica- 
tive calibration changes of 1.02, 1.05 and 0.94 at 250, 350 and 
500 //m respectively relative to the previous Ceres flux calibra- 
tion (Griffin et al. 2010, Swinyard et al. 2010). Jumps in ther- 
mometer timelines were identified and corrected using the mean 
levels on either side of each jump. Thermal drift modelling was 
achieved through low-pass filtering. Bolometer timelines were 
filtered using high-pass 4 mHz filtering to remove correlated 1 // 



noise. Maps were made using the naive map maker in HIPE and 
filtered optimally for point sources. 

2.2. Quasar compilations 

The survey areas of H-ATLAS were chosen for their supplemen- 
tary multi-wavelength coverage. The H-ATLAS science demon- 
stration field is covered by the SDSS with 240 QSOs in its spec- 
troscopic QSO catalogue in the field. There are also 171 QSOs 
from the 2dF-SDSS LRG and QSO (2SLAQ) survey. 

We supplement our H-ATLAS QSOs with the compilation 
from Serjeant & Hatziminaoglou (2009), including the authors' 
IRAS ADDSCAN photometry of Palomar-Green QSOs and 
5/?/fzer photometry of QSOs from the Spitzer SWIRE legacy sur- 
vey (Lonsdale et al. 2003, 2004), together with QSOs observed 
at 850 yum or 1200yum (Omont et al. 1996, Omont et al. 2001, 
Cai-ifli et al. 2001, Isaak et al. 2002, Omont et al. 2003, Priddey et 
al. 2003a, 2003b, Robson et al. 2004, Wang et al. 2007). Figure 
[1] shows the distribution of our QSO compilation in the opti- 
cal luminosity-redshift plane. Note that by combining several 
flux-limited samples we mitigate the effects of Malmquist bias 
(Teerikorpi 1984), i.e. that in a single flux-limited sample there 
is a strong degeneracy between luminosity and redshift, making 
it impossible to decouple evolution effects from luminosity de- 
pendence. For the time being we use the working assumption 
that any subtle biases in e.g. the Palomar-Green sample are not 
strongly correlated with far-infrared luminosity; in time the full 
H-ATLAS QSO sample will allow us to make use of a more uni- 
form QSO selection. 



3. l\/letliods 

We adopt the methodology of Serjeant et al. (2004, 2008, 2009) 
for comparing QSO fluxes with the flux distributions of the H- 
ATLAS maps as a whole. We use unweighted means to estimate 
H-ATLAS QSO fluxes and use the Kolmogorov-Smirnoff statis- 
tic to compare flux distributions. The advantage of this statistic 
is that it contains its own control test, i.e. there is no need for 
performing stacks on randomised QSO positions. 

While most QSOs can only be detected in the far-infrared 
and submm through stacking analyses, there are a few QSOs 
with direct detections. How should one combine these? For ex- 
ample, if one has a single non-detection of + 4 mJy and single 
detectionof 1 00+4 mJy, what can one say about the average flux 
fi of this population? Clearly the answer is not/i - 50±2 mJy. We 
have adopted the methodology of Serjeant and Hatziminaoglou 
(2009) and treated flux measurements of individual QSOs as 
attempts to determine the population mean. The dispersion in 
the population is an error term on this measurement, which one 
would add in quadrature to the noise on any individual measure- 
ment. We determine the population dispersion from our data and 
made a simultaneous maximum-likelihood fit to the mean ju and 
dispersion cr. There is no covariance between these parameters 
because of the independence of signal and noise. Where there are 
fewer than 3 objects being considered we set cr = 0.84yU. More 
details can be found in Serjeant & Hatziminaoglou (2009). 

We used this method to combine direct detections with non- 
detections. We also combine the Serjeant & Hatziminaoglou 
(2009) sample with our H-ATLAS data. We use the point- 
source-convolved H-ATLAS maps (Pascale et al. 2010). 
Measurements taken off" these maps report the point source flux 
plus a background flux contribution from other galaxies. To set 
the latter to zero we convolve our maps with a further kernel K 
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with a zero total, i.e. J K(r)27Trdr - 0. K was set to a constant 
negative value at radii r of 2.5 to 6 times the point spread func- 
tion FWHM, with a unit S function in its central pixel, and zero 
everywhere else. After convolution with K, every galaxy makes 
an exactly zero net contribution to the map. If there is a non-zero 
angular cross-correlation function from galaxies with a physi- 
cal association with the QSO, then they may still contribute to 
the far-infrared and submm fluxes. One approach is to use con- 
straints on the angular correlation function of these galaxies to 
place a bound on this contribution (e.g. Serjeant et al. 2008). 
However, in this case we are studying the assembly of the stel- 
lar mass and the companion galaxies may in time be accreted by 
the QSO host galaxy. We have therefore chosen to associate all 
the star formation inferred from the flux in the far-infrared beam 
with the QSO host galaxy. 

4. Results 

4.1. Stacking analysis resuits 

We compared the flux measurements at the positions of our 
H-ATLAS QSOs with the distribution of flux measurements 
throughout the H-ATLAS maps. The Kolmogorov-Smirnoff' 
statistic rejects the null hypothesis that these are drawn from the 
same distribution at significance levels of 3 x 10""'^, 7 x 10"^° 
and 7 x lO"** at 250, 350 and 500 ;um respectively, equivalent to 
12.0cr, 9.1cr and 5.7cr. This is not due simply to the presence 
of bright sources, since removing QSOs with fluxes > lOOmJy 
and comparing with regions of the map with flux measurements 
< 100 mJy still yields significance levels of 11.9cr, 9.0cr and 
5.6cr. The average fluxes of the QSOs in the H-ATLAS sci- 
ence demonstration field are S250fim - 9.41 + 0.88 mJy (llcr 
detection), ^sgo^m - 7.68 + 0.87 mJy (8.9cr detection) and 
Ssoofim = 5.14 + 0.92 mJy (5.6cr detection). 

When calculating lOO/^m rest-frame luminosities for H- 
ATLAS QSOs, we use the closest SPIRE filter to 100(1 + 
z)i-im. Table [1] lists our estimates of the mean lOOjum rest- 
frame luminosities of the whole QSO compilation in redshift 
and optical luminosity bins, using the methodology of Sect. 
[3] The results are well-fit (xl - 0.93) by the expression 
logio(vLv(lOO;um)/lO'2L0) = a(z)[lAB +/3(z)], where a(z) = 
(0.0371 + 0.0048) X min(z,4) - 0.235 + 0.018, y6(z) = (-1.19 + 
0.30) X min(z,4) + 27.42 + 0.37, and min(z,4) = z at z < 4 
and 4 otherwise. The slope of the luminosity-luminosity cor- 
relation is shallower at high redshift, as found by Serjeant & 
Hatziminaoglou 2009 (see also Mullaney et al. 2010). 

QSO number densities 'S>{Iab,z) at these absolute magni- 
tudes and redshifts are already well-determined. We adopted the 
luminosity functions of Croom et al . 2004 at z < 2 . 1 and the pure 
density evolution model of Meiksin 2005 at z > 3, with an opti- 
cal spectral index of -0.5 to transform to /-magnitudes. Between 
z = 2.1 and z = 3 we interpolate in log <l> between these mod- 
els at fixed optical luminosities. We assume the far-infrared is 
dominated by giant molecular clouds so vLy(lOOjum)/lO'^L0 = 
265Mo/year as appropriate for our assumed SED and a Salpeter 
initial mass function from 0. 1 to 100 Mq (Kennicutt 1998). Table 
[1] also lists the comoving star formation densities inferred from 
vLv(lOOyum) X <I>(Iab,z). 

Figure |2] plots the data in table [T] interpolating between the 
upper and lower bounds of the comoving star formation densi- 
ties at the mid-points of each redshift bin. The evolution in the 
comoving star formation density of -22 > /ab > -24 QSOs 
from the 0.05 < z < 0.5 bin to the 1 < z < 2 bin scales as 
(1 -I- z)^"'"'"'^, but the evolution at higher luminosities is much 
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Table 1. Average stacked QSO 100 /im rest-frame luminosities 
for bins in redshift and /-band absolute magnitude. An M82 SED 
has been assumed, though the results are only very weakly de- 
pendent on the assumed SED. The first column gives the number 
of QSOs in the bin in question. The final column gives the co- 
moving volume-averaged star formation densities of QSO host 
galaxies in each bin. 



faster: at -24 > /ab > -26 the variation is (1 -H z)^'^='""', while 
at -26 > /ab > -28 it reaches an astonishing (1 + z)'"''' over 
this redshift range (mostly but not entirely due to the luminos- 
ity function). There are no QSOs at /ab < -28 at z < 0.5 but 
between the 0.5 < z < 1 and 1 < z < 2 bins, the evolution of 
the volume-averaged star formation in /ab < -28 QSOs scales 
as (1 -H z)^^^. An important caveat is that we are only address- 
ing the optically-defined QSO population, so we are necessarily 
missing the type-2 QSOs and Compton-thick objects (see e.g. 
the discussion in Almaini, Lawrence & Boyle 1999). A simi- 
lar downsizing effect is seen in the QSO soft X-ray luminosity 
function (e.g. Hasinger et al. 2005). The type-2 QSO fraction in- 
creases with redshift (e.g. Hasinger 2008) which would increase 
our inferred evolution rates. 



5. Discussion and conclusions 

Could the co-evolution of the total cosmic star formation his- 
tory and total black hole accretion be explained by both being 
simultaneously driven by major galaxy-galaxy mergers (e.g. di 
Matteo et al. 2005)? Most star formation at z < 1 does not seem 
to be triggered by major mergers (e.g. Bell et al. 2005), but QSO 
host galaxies often have signs of significant disturbance. This 
suggests a more complex picture. An alternative possibility pro- 
posed by Zheng et al. (2009) is that star formation occurs mainly 
in disks, while black hole growth occurs separately after ma- 
jor merger events. Any apparent co-evolution in this model is 
present because the black hole accretion depends on the potential 
well of the spheroid, which in turn grows through the addition 
of disrupted disk stars. 

Our observation of an unprecedented strong evolution in the 
star formation rates of the brightest QSOs has no obvious ana- 
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Fig. 2. Cosmic star formation history of QSO host galaxies inferred from 100 //m rest-frame luminosities, for QSOs with -22 > 
/ab > -24 (red), -24 > /ab > -26 (orange), -26 > /ab > -28 (green) and /ab < -28 (blue). The 2 < z < 4, -22 > /ab > -24 data 
point has too high a noise level to be usefully constraining and has been omitted for clarity. For comparison, the z = total galaxy 
star formation rate is (2.9 + 0.7) x IO^^Mq yr"' Mpc"^ (e.g. Serjeant et al. 2002). 



logue in evolving galaxy populations. This is difficult to accom- 
modate in models with quasi-simultaneous starburst and AGN 
events sharing a common trigger. Instead, it suggests star for- 
mation in bright QSO hosts is not representative of stellar mass 
assembly in general. We find QSO hosts comprise i 1% of the 
total volume-averaged star formation density of the Universe at 
all redshifts (e.g. Michalowski et al. 2010). 

What could drive the extremely strong downsizing evolution 
in bright QSOs? If major mergers drive QSO activity, the declin- 
ing major merger rate at a given halo mass could account for the 
decreasing QSO number density. However, this is not enough to 
account for our observations, since we also find that the relation- 
ship between optical luminosities (our proxy for black hole ac- 
cretion rate) and far-infrared luminosities (our proxy for star for- 
mation rate) changes with redshift, in the sense that bright QSOs 
at a fixed optical luminosity were more far-infrared-luminous at 
high redshifts. One possibility is that less gas is available at low 
redshift for star formation and black hole accretion, due to con- 
sumption in star formation or ejection through feedback mech- 
anisms. The strong decline in the comoving star formation den- 
sities of bright QSOs could therefore be due to a double depen- 
dence on the availability of gas combined with a declining major 
merger rate. 
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